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Optimization Methods for Designing Sequences 
with Low Autocorrelation Sidelobes 
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Abstract —Unimodular sequences with low autocorrelations are 
desired in many applications, especially in the area of radar 
and code-division multiple access (CDMA). In this paper, we 
propose a new algorithm to design unimodular sequences with 
low integrated sidelobe level (ISL), which is a widely used 
measure of the goodness of a seqnence’s correlation property. 
The algorithm falls into the general framework of majorization- 
minimization (MM) algorithms and thns shares the monotonic 
property of such algorithms. In addition, the algorithm can be 
implemented via fast Fourier transform (FFT) operations and 
thus is computationally efficient. Fnrthermore, after some mod¬ 
ifications the algorithm can be adapted to incorporate spectral 
constraints, which makes the design more flexible. Numerical 
experiments show that the proposed algorithms outperform 
existing algorithms in terms of both the qnality of designed 
seqnences and the compntational complexity. 

Index Terms —Unimodular sequences, autocorrelation, inte¬ 
grated sidelobe level, majorlzatlon-minimization. 

I. Introduction 

Since the 1950s, digital communications engineers have 
sought to identify sequences whose autocorrelation sidelobes 
are collectively as low as possible according to some suitable 
measure of “goodness”. Applications of sequences with low 
autocorrelation sidelobes range from synchronization to code¬ 
division multiple access (CDMA) systems and also radar sys¬ 
tems HI, 0. Low autocorrelation can improve the detection 
performance of weak targets El in range compression radar 
and it is also desired for synchronization purposes in CDMA 
systems. Additionally, due to the limitations of sequence 
generation hardware components (such as the maximum sig¬ 
nal amplitude clip of analog-to-digital converters and power 
amplifiers), it is usually more desirable to transmit unimodular 
(i.e., constant modulus) sequences to maximize the transmitted 
power available in the system a. 

The focus of this paper is on the design of unimodular 
sequences with low (aperiodic or periodic) autocorrelation 
sidelobes. Let {xn]n=i denote the complex unimodular (with¬ 
out loss of generality, we will assume the modulus to be 1) 
sequence to be designed. The aperiodic (r^) and periodic (f^) 
autocorrelations of {xn}n=i ^re defined as 

N-k 

X =r%,k = Q,...,N -I, 

n—1 
N 

X! (n+k) (mod N) ~ ^-fe) ^ = 0, . . . , TV — 1 . 

n—1 
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Then the integrated sidelobe level (ISL) of the aperiodic 
autocorrelations is defined as 

Af-l 

ISL=^|rfe|2, (1) 


which will be used to measure the goodness of the aperiodic 
autocorrelation property of a sequence. The ISL metric for the 
periodic autocorrelations can be defined similarly. Note that 
the ISL metric is highly related to another important goodness 
measure; the merit factor (MF), which is defined as the ratio 
of the central lobe energy to the total energy of all other lobes 
by Golay in 1972 Q, i.e.. 


MF = 




2ISL’ 


( 2 ) 


Owing to the practical importance and widespread appli¬ 
cations of sequences with good autocorrelation properties, in 
particular with low ISL values or large MF values, a lot of ef¬ 
fort has been devoted to identifying these sequences via either 
analytical construction methods or computational approaches 
in the literature. The studies focused more on the design of 
binary sequences at the early stage, and extended to polyphase 
sequences later; see 0-iini for the case of binary sequences 
and mi-iiisi for the case of polyphase or unimodular se¬ 
quences. Some sequences with good autocorrelation properties 
that can be constructed in closed-form have been proposed in 
the literature, such as the Frank sequence lfT2l and the Golomb 
sequence E]. Computational approaches, such as exhaustive 
search IS), evolutionary algorithms 0, heuristic search ani 
and stochastic optimization methods m, na, have also been 
suggested. However, these methods are generally not capable 
of designing long sequences, say N ~ 10^ or larger, due to the 
increasing computational complexity. Recently, an algorithm 
named CAN (cyclic algorithm new) was proposed in lEl 
that can be used to produce unimodular sequences of length 
N oo 10® or even larger with low aperiodic ISL and it has 
been adapted to generate sequences with impulse-like periodic 
autocorrelations im and incorporate spectral constraints ll20l . 
But rather than the original ISL metric, the CAN algorithm 
tries to minimize another simpler criterion which is stated to 
be “almost equivalent” to the ISL metric. 

In this paper, we develop an efficient algorithm that directly 
minimizes the ISL metric (of both aperiodic and periodic auto¬ 
correlations) monotonically, which we call MISL (monotonic 
minimizer for integrated sidelobe level). The MISL algorithm 
is derived via applying the general majorization-minimization 
(MM) method to the ISL minimization problem twice and 
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admits a simple closed-form solution at every iteration. The 
convergence of the algorithm to a stationary point is proved. 
Similar to the CAN algorithm, the proposed MISL algorithm 
can also be implemented via fast Fourier transform (FFT) 
operations and is thus very efficient. But due to the nature of 
double majorization, the MISL algorithm may converge very 
slow especially for large N, and we propose to apply some 
acceleration schemes to fix this issue. A MISL-like algorithm, 
which shares the simplicity and efficiency of MISL, has also 
been developed for minimizing the ISL metric when there are 
additional spectral constraints. 

The remaining sections of the paper are organized as 
follows. In Section II, the problem formulations are presented 
and the CAN algorithm is reviewed. In Section III, we first 
give a brief review of the MM framework and then the MISL 
algorithm is derived, followed by the convergence analysis. In 
Section IV, some acceleration schemes for MISL are presented 
and then the spectral-MISL algorithm is developed in Section 
V. Finally, Section VI presents some numerical examples and 
the conclusions are given in Section VII. 

Notation: Boldface upper case letters denote matrices, bold¬ 
face lower case letters denote column vectors, and italics 
denote scalars. R and C denote the real field and the complex 
field, respectively. Re(-) and Im(-) denote the real and imag¬ 
inary part respectively, arg(-) denotes the phase of a complex 
number. The superscripts (•)^, (•)* and (•)^ denote transpose, 
complex conjugate, and conjugate transpose, respectively. Xij 
denotes the (/-th, y-th) element of matrix X and Xi denotes the 
/-th element of vector x. Tr(-) denotes the trace of a matrix. 
diag(X) is a column vector consisting of all the diagonal 
elements of X. Diag(x) is a diagonal matrix formed with x 
as its principal diagonal. vec(X) is a column vector consisting 
of all the columns of X stacked. I„ denotes an n x n identity 
matrix. 


where ujp = ^{p — 1), P = I,'" - , ‘2N. Thus, the sequence 
design problem Q becomes 


minimize 

Xn 

2N 

E 

N 

1- 

1 

CM 


p^i 

n—1 


subject to 

\Xn\ 

= 1, n= 1,...,TV. 


Let us define 

X = [xi, 

■■ ,Xn]'^ 



n 2 


(5) 


(6) 


a 


ip 




then the problem (|5]l can be rewritten as 

2N 

minimize [a^xx^ap — 7V]" 


p=i 


( 8 ) 


subject to |a;n,| = 1, n = 1,..., iV. 
Expanding the square in the objective yields 


minimize 


2N 

E 

p^i 


a^xx-^aj 




■p 




subject to \xn\ = 1, n = 1,... ,N, 


(9) 


where the second term in the objective can be shown to be 
a constant (using ParsevaTs theorem), i.e., = 

2A^||x||2 = 27V^. Thus by ignoring the constant terms, the 
problem can be simplified as 

2N 

minimize > (a^xx-^ap) 

X 2 —/ ' F F/ (iU) 

p=l 

subject to \xn\ = n = 1 ,..., N. 

The problem is hard to tackle, due to the quartic objective 
function and also the nonconvex unit-modulus constraints. 


B. The Periodic Autocorrelation 


IT Problem Formulation and Existing Methods 


The problem of interest is the design of a complex uni- 
modular sequence {xn}n=i that minimizes the ISL metric (of 
either aperiodic or periodic autocorrelations), i.e.. 


minimize ISL 

subject to |a;n,| = 1, n = 1,..., iV. 


(3) 


In the following, we first reformulate the problem Q for the 
case of aperiodic and periodic autocorrelations respectively 
and then briefly review two corresponding algorithms pro¬ 
posed in the literature. 


A. The Aperiodic Autocorrelation 


Similarly, in the periodic case, the ISL metric ISL = 
SE/ also be expressed in the frequency domain 

as 113 

2 


N 


ISL = - y 
N ^ 




N 

2 


-N 

n—1 



( 11 ) 


where Cdp = ^{p — 1), p = 1, • • • , A^. By defining 

ap = (12) 

the ISL minimization problem Q in this case can be rewritten 


minimize 


N 

y [afxx"ap-7V]' (13) 

p=i 

subject to |a;„| = 1, n = 1,..., iV. 


It has been shown in iflTl that the ISL metric of the aperiodic 
autocorrelations can be expressed in the frequency domain as 


2N 

'"■'“iivi: 




N 

2 


-N 

n—1 



Eollowing similar steps as in the previous subsection, the 
problem can be further simplified as 

N 

minimize (a^xx-^ap) 

p=i 

subject to |a;n,| = 1, n = 1,..., TV. 


(4) 


(14) 
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We note that although problem (fT4l i looks similar to the 
problem ( fTOl i in the aperiodic case, it is usually a simpler task 
to find sequences with good periodic autocorrelation proper¬ 
ties. In particular, in the periodic case unimodular sequences 
with zero ISL exist for any length N mi, however it is 
not possible to construct unimodular sequences with exact 
impulse-like aperiodic autocorrelations. 


C. Existing Methods 

An algorithm named CAN was proposed in ifTTl for de¬ 
signing unimodular sequences with low ISL of aperiodic 
autocorrelations. But instead of minimizing the ISL metric 
directly, i.e., solving Q, the authors proposed to solve the 
following simpler problem 


minimize 
subject to 


N 


2N 

E 

p—l n—l 

\xn\ = l,n = l,...,N, 




(15) 


whose objective is quadratic in {xn}- Let A be the following 
N X 2N matrix; 


A — [ai,..., a2Ar], 


(16) 


then the objective function of problem (fTST i can be rewritten 
in the following more compact form: 


A^x — VNv 


(17) 


where v = i The CAN algorithm then 

minimizes ( fTTl ) by alternating between x and v. For a given 
X, the minimization of ( fTTl l with respect to {V'p} yields 


V'p = arg(/p), p = 1,.. .,2A, (18) 

where f = A^x. Similarly, for a given v, the minimizing x 
is given by 

Xn = n = 1,..., N, (19) 


where g = Av. The CAN algorithm is summarized in 
Algorithm [T] for ease of comparison later. 


Algorithm 1 The CAN algorithm proposed in ifTTl . 
Require: sequence length N 
1: Set k = 0, initialize x(°l. 

2: repeat 

3: f = A^xl'^l 

4: Up = e^“®(/p),p = l,...,2A 

5: g = Av 

6: n= 

7: A: ^ fc -I- 1 

8: until convergence 


It was mentioned in IfTTl that the problem ( fTSl ) is “almost 
equivalent” to the original problem Q. But as the authors 
also pointed out, the two problems are not exactly equivalent 


and they may have different solutions. It is easy to see that 
minimizing with respect to {^’p} in ( flST l leads to 


2N ■ 

N 

" 

minimize > 

Xn * ^ 


- y/N 

. 

n—l 

_ 


subject to \xn\ = 1, n = 1,..., A, 


which is different from the original problem (|5]l. So the 
point the CAN algorithm converges to is not necessary a 
local minimum (or even a stationary point) of the original 
ISL metric minimization problem. Although CAN does not 
directly minimize the ISL metric, in IfTTl the authors showed 
numerically that CAN could generate sequences with good 
correlation properties. Moreover, CAN can be implemented 
via FFT, thus can be used to design very long sequences. 

In the periodic case, an algorithm called PeCAN was used 
in in to construct unimodular sequences with low ISL. It 
can be viewed as an adaptation of the CAN algorithm for the 
periodic case and similarly instead of minimizing the original 
ISL metric (fTTT) . the PeCAN algorithm considered the simpler 
criterion 


E 


N 


n—l 


2 


( 21 ) 


where {'0p} are auxiliary variables as in ( fTSl . It was shown 
numerically in lfT9l that PeCAN can generate almost perfect 
(with zero ISL) sequences from random initializations. 

Although the criteria considered by CAN and PeCAN 
are somehow related to the original ISL metric and the 
performance of the two algorithms is acceptable in some 
situations, one may expect that directly solving the original 
ISL minimization problem can probably lead to better perfor¬ 
mance. In the next section, we will develop an algorithm that 
directly minimizes the ISL metric, while at the same time is 
computationally as efficient as the CAN algorithm. 


III. ISL Minimization via 
Majorization-Minimization 

In this section, we first introduce the general majorization- 
minimization (MM) method briefly and then develop a simple 
algorithm for the ISL minimization problem based on the 
general MM scheme. In the derivation, we will focus on 
the aperiodic case, i.e., the problem (fTOl i. but the resulting 
algorithm can be readily adapted for the periodic case. 


A. The MM Method 

The MM method refers to the majorization-minimization 
method, which is an approach to solve optimization problems 
that are too difficult to solve directly. The principle behind 
the MM method is to transform a difficult problem into a 
series of simple problems. Interested readers may refer to lITTI 
and references therein for more details (recent generalizations 
include 1221, CD). 

Suppose we want to minimize /(x) over X G C". 
Instead of minimizing the cost function /(x) directly, the 
MM approach optimizes a sequence of approximate objective 
functions that majorize /(x). More specifically, starting from a 
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feasible point the algorithm produces a sequence 
according to the following update rule 


necessary. Let us define X = xx^ and Ap = a^a^, then the 
problem (fTOl i can be rewritten as 


x('=+i) e argmin w(x,x('')), (22) 

xex 

where x^^^ is the point generated by the algorithm at iteration 
k, and u (x, x^^^) is the majorization function of /(x) at 
x(^). Formally, the function u(x,xW) is said to majorize the 
function /(x) at the point provided 

u(x,x('')) > /(x), \/xGX, (23) 

u(x('=\xW) = /(xW). (24) 

In other words, function m(x, x^^^) is an upper bound of /(x) 
over X and coincides with /(x) at 

To summarize, to minimize /(x) over X G C", the main 
steps of the majorization-minimization scheme are 

1) Find a feasible point x^^^ and set k = 0. 

2) Construct a majorization function tt(x, x^^^) of /(x) at 
x(fc)^ 

3) Let x(^+^^ G argmin u(x, 

xex 

4) If some convergence criterion is met, exit; otherwise, set 
k = k + 1 and go to step (2). 

It is easy to show that with this scheme, the objective value 
is decreased monotonically at every iteration, i.e., 

/(x(fc+i)) < u(x('=+i),xW) < u(xW,xW) = /(xW). 

(25) 

The first inequality and the third equality follow from the the 
properties of the majorization function, namely (l23t and (l24l i 
respectively and the second inequality follows from (l22l i. The 
monotonicity makes MM algorithms very stable in practice. 


B. MISL 

To solve the aperiodic problem (fTOl i via majorization- 
minimization, the key step is to hnd a majorization function 
of the objective such that the majorized problem is easy to 
solve. For that purpose we first present a simple result that 
will be useful when constructing the majorization function. 

Lemma 1. Let h be an n x n Hermitian matrix and M be 
another nx n Hermitian matrix such that M ^ L. Then for 
any point xq G C", the quadratic function x^Lx is majorized 
by x^Mx + 2Re (x^(L — M)xo) + x^(M — L)xo at xq. 

Proof: It is easy to check that the two functions are equal 
at point Xg. Since M ^ L, we further have 

x^Lx 

= x^Lxo-|-2Re ((x - xo)-^Lxo) -l-(x - xo)-^L(x - xg) 

< x^Lxg-|-2Re ((x - Xgj'^Lxg) +(x - Xo)'^M(x - Xg) 
= x^Mx+2Re (x^(L — Mjxg) -f x^(M — Ljxg 

for any x G C”. The proof is complete. ■ 

The objective of the problem ( fTOl i is quartic with respect to 
X. To hnd a majorization function, some reformulations are 


minimize 

x,X 

subject to 


2N 

^Tr(XA,)2 

p=i 

X = xx^ 

\xn\ = 1, n = 1,... ,7V. 


(26) 


Since Tr(XAp) = vec(X)'^vec(Ap), problem (l26l l is just 


minimize vec(X)^€>vec(X) 

subject to X = xx^ 

|a;„| = 1, n = 1,..., A, 


(27) 


where $ = vec(Ap)vec(Ap)^. We can see that the 

objective of dZTl l is a quadratic function in X now. Given 
X(fe) = at iteration k, by choosing M = 

Amax(^)I in Lemma [T] we know that the objective of (l27T i 
is majorized by the following function at 

wi(X,X«) 

=Amax(^)vec(X)^vec(X) 

+ 2Re (vec(X)^($ - A^ax(^)I)vec(X('=))) 

+ vec(X('=) (A„,ax(^)I - $)vec(X('=)) 


where Aniax(^) is the maximum eigenvalue of $ and can be 
shown to be 2A^ (see Appendix lAli. Since vec(X)^vec(X) = 
(x^x)^ = A^, the hrst term of (l28T l is just a constant. After 
ignoring the constant terms in (l28T l. the majorized problem of 
(l27l i is given by 

minimize Re (vec(X)^ (# - 2A2 i) vec(X('=))) 

subject to X = xx^ (29) 

\xn\ = 1, n= 1,...,A, 


which can be rewritten as 


minimize 

x,X 

subject to 


2N 

^Tr(X('=)Ap)Tr(ApX) - 2A2Tr(X('=)X) 

p=i 

X = xx^ 

\xn\ = 1, n = 1,..., A. 

(30) 


Remark 2. To solve the problem (l30l l. one possible approach 
is to apply an SDP relaxation, yielding 


2N 


minimize y]]Tr(X('=)Ap)Tr(ApX) - 2A2Tr(XWX) 


subject to X ^ 0 


diag(X) = 1. 


(31) 

The above SDP is empirically observed to always give a rank 
one solution, but proving the tightness of the SDR is out of the 
scope of this paper. In addition, the computational complexity 
of solving an SDP is high and it needs to be solved at every 
iteration, which makes it not amenable for large A. 


In the following, we propose to apply majorization- 
minimization to the problem (l30l l again, which will yield a 
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simple algorithm. To apply the second majorization, we hrst 
rewrite the problem ( l30l l as 


2N 


minimize P la^xP — 2iV^ Ix^x^^^ 1^ 

X ' P I I P I I I 

p=i 

subject to |xn| = 1, n = 1,..., TV. 

It can be written more compactly as 

(ADiag(p('=))A'f^ - 2A2 x('=) (x^'^)) 


(32) 


H 


minimize x 

X 

subject to |x„| = 1, n = 1,..., A, 

(33) 

where A = [ai,..., a 2 Ar] and (Ij^ denotes 

the element-wise absolute-squared value). We can clearly see 
that the objective function in (l3Tt is quadratic in x and by 
choosing M = pmlxAA^ in Lemma [1] it is majorized by the 
following function at x^^^; 


U2(x,xW) 


=nl5^Lx^AA^x + 2Re 


(x«(A-2A2xW(xW)")xW) 


(k) 


+ (x('=))^(2A2x('=)(xW)^ - A)x' 

J34) 

where pmlx = maxpjpp^^ : p = 1,...,2N} and A = 


(^Diag(p('=)) 




■ Pmaxlj A^ . Since AA^ = 2AI, the first 
term of ( l34l i is a constant. By ignoring the constant terms in 
(l34l i. we have the majorized problem of (l33T i: 


minimize Re 


(x^ (A-2A2 x«(x«)^) xW) 


subject to |x„| = 1, n = 1,..., A, 
which can be rewritten as 

minimize ||x —y ||2 


subject to |x„| = 1, n = 1,..., A, 


where 


y = - (a - 2A2x('=) 


Ak) 


= -A 


(Diag(p('=)) -p«xl - 


(35) 


(36) 


(37) 


It is easy to see that the problem ( l36l l has a closed-form 
solution, which is given by 

x„ = e^"s(^"\n = l,...,A. (38) 

Now we can summarize the overall algorithm and it is given 
in Algorithm |2] Note that although in the derivation we have 
applied the majorization-minimization scheme twice, it can be 
viewed as directly majorizing the objective of the problem (fTOl i 
at x^^) over the constraint set by the following function; 


u(x,xW) 


2N 


=2u2(x, -I- AN* - a^x 


(k) 


p=i 


=4Re (x^A (Diag(p('=)) - - A^l) A^x^''^) 

2N . 

+ 8A2(pW, + A2) -3y]|a^x«| 


p=i 


and the algorithm preserves the monotonicity of the general 
majorization-minimization method. Thus, we call the algo¬ 
rithm MISL (Monotonic minimizer for Integrated Sidelobe 
Level). 

We also note that the above derivations can be carried out 
similarly for the periodic problem (fl4l l and it turns out that 
the MISL algorithm can be easily applied for the periodic 
case after replacing A in (fTbl l and Algorithm |2] with A = 
[ai,..., ajv], where {hp}^]^ are defined in (fT^. 


Algorithm 2 MISL - Monotonic minimizer for Integrated 
Sidelobe Level._ 

Require: sequence length A 
1: Set k = 0, initialize x*^°\ 

2: repeat 

4: Pmlx = maxpjp^''^ : p = 1,..., 2A} 

5: y = - A (Diag(p('=)) - pLIxI - A^l) A^x^'^) 

6: ri= 1,..., A 

7: k-l^ k+1 

8: until convergence 


C. Convergence Analysis 

The MISL algorithm given in Algorithm |2] is based on the 
general maximization-minimization scheme, thus according to 
subsection IIII-AI we know that the sequence of objective 
values evaluated at generated by the algorithm is 

nonincreasing. And it is easy to see that the objective value is 
bounded below by 0, thus the sequence of objective values is 
guaranteed to converge to a hnite value. 

In this section, we will further analyze the convergence 
property of the sequence {x^^^} generated by the MISL 
algorithm and show the convergence to a stationary point. 

We hrst introduce a hrst-order optimality condition for 
minimizing a smooth function over an arbitrary constraint set, 
which follows from ll24l . 


Proposition 3. Let f : R" H be a smooth function, and 
let X* be a local minimum of f over a subset X o/R". Then 

V/(x*)^y > 0, Vy G T;t(x*), (40) 

where Tx{'x*) denotes the tangent cone of X at x*. 


A point X G A satisfying the hrst-order optimality condition 
(l40l i will be referred as a stationary point. 

Next we exploit a nice property shared by the minimization 
problems under consideration. 


Lemma 4. Let g : ^ R he a real valued function with 

complex variables. Then the problem 


minimize 
xec" 
subject to 

is equivalent to 


minimize 
xec" 
subject to 


5(x) 

\xn\ = 1, n = 1,..., A 

p(x) + ax^x 

\xn\ = 1, n= 1,...,A 


(41) 


4 


( 42 ) 
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for any finite value a, in the sense that they share the same 
optimal solutions. 

Proof: Since the additional term ax^x is just the constant 
aN over the constraint set {x € C^| \xn\ = l,n = 
1,...,7V}, it is obvious that any solution of problem (ITTI l 
is also a solution of problem (l42l i and vice versa. ■ 

To facilitate the analysis, we further note that upon defining 


where 

2N 

p=i 

(49) 

Thus, x(°°) = [Re(x(°°^)^, Im(x(°“^)^]'^ is a global mini- 
mizer of (l48T l. Then as a necessary condition, we have 


Ap — 


X = [Re(x)'^, Im(x)'^]'^ 
Re(Ap) —Im(Ap) 


Im(Ap) Re(Ap) 


(43) 

(44) 


it is straightforward to show that the quartic complex mini¬ 
mization problem ( fTOb is equivalent to the following real one; 


2N , 

linimize ^x^ApX^ 


p=i 

subject to x: 


(45) 


' xl^N = l,n = l,...,N. 


We are now ready to state the convergence properties of 
MISL. 


Theorem 5. Let be the sequence generated by the 

MISL algorithm in Algorithm |2] Then every limit point of the 
sequence is a stationary point of the problem (Hoj. 

Proof: Denote the objective functions of the problem (fTOl i 
and its real equivalent (l45l l by /(x) and /(x), respectively. De¬ 
note the constraint sets of the problem (fTOl) and (l45l l by C and 
C, respectively, i.e., C = {x S C^| |x„| = 1, n = 1,..., N} 
and C = {x G = 1, n = 1,..., N}. Then 

from the derivation of MISL, we know that /(x) is majorized 
by the function rt(x, in ( |39] | at x^^^ over C. According 
to the general MM scheme described in subsection IIII-AI we 
have 


/•(x(fc+i)) < u(x('=+i),xW) < m(xW,xW) = /(xW), 


which means {/(x^*^))} is a nonincreasing sequence. 

Assume that there exists a converging subsequence 
x(“), then 

+ + = /(x('=^ + l)) < /(x('=^+l)) 

< u(x('=^+i),x('=^4) < m(x,x('=^)),Vx G C. 
Letting j —>■ -|-c», we obtain 

x^°°^) < m(x, x^°“^), Vx G C, (46) 


i.e., x(°“) is a global minimizer of u(x, x(°“)) over C. Accord¬ 
ing to Lemma |4] x(°°) is also a global minimizer of 


minimize u(x, x^°°^)-|-4A^x^x 

X 

subject to X G C. 


(47) 


With the definitions of x and Ap given in (l43T l and (l44l i and 
by ignoring some constant terms in m(x, x(°“)), it is easy to 
show that (l47l i is equivalent to 


2N 

minimize 2p^^iy^x^ApX-b 4x^b-f 4A^x^x 
subject to X G C, 


(48) 


Vu(5c(“))^y>0, VyGrc(k(“)), (50) 

where u{x) denotes the objective function of (l48T l. It is easy 
to check that 

2N 

= Vu(k(°°)) = 4^ ((x(°°))^Apx(°°)) Apx(“). 


Thus we have 

V/(x(-))ry > 0, Vy G Tc(x(°“)), (51) 

implying that is a stationary point of the problem (l45T l. 
Due to the equivalence of problem (l45l l and (fTOl i. the proof is 
complete. ■ 


D. Computational Complexity of MISL 

From Algorithm 121 we can see that the per iteration compu¬ 
tational complexity of MISL is dominated by two matrix vec¬ 
tor multiplications involving A, which can be easily computed 
via FFT and IFFT operations. More specifically, let F be the 
2N X 2N DFT matrix with F^ n = , 0 < to, n < 2N, 

it is easy to see that the N x 2N matrix A is composed of 
the first N rows of F^. Given a vector x G C^, the product 
A^x can be computed via the FFT of the zero padded version 
of X, i.e., Flx-^jO^]-^. Similarly, given a vector z G 
the product Az consists of the first N elements of the inverse 
FFT of z (i.e., F^z). In the periodic case, since the matrix 
A^ itself is the N x N DFT matrix, multiplying a vector by 
and A is just the FFT and IFFT of the vector. Thus, the 
MISL algorithm is computationally very efficient and can be 
used for the design of very long sequences, say N ^ 10®. 


IV. Acceleration Schemes 

As described in the previous section, the derivation of 
MISL is based on majorization-minimization principle, and 
the nature of the majorization functions will dictate the con¬ 
vergence speed of the algorithm. From numerical simulations, 
we noticed that for large N, the convergence of MISL is very 
slow, which may be due to the double majorization scheme 
that we carried out in its derivation. One option to fix this 
issue is to employ some acceleration schemes to accelerate 
the convergence of MISL. In the following, we will consider 
two such schemes (for the aperiodic case) and the resulting 
algorithms can be adapted for the periodic case by using A 
instead of A. 
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A. Acceleration via Fixed Point Theory 

In this subsection, we describe an acceleration scheme that 
can be applied to accelerate the MISL algorithm proposed 
in this paper. It is the so called squared iterative method 
(SQUAREM), which was originally proposed in ll25l to ac¬ 
celerate any EM algorithms. SQUAREM follows the idea 
of the Cauchy-Barzilai-Borwein (CBB) method Il26ll . which 
combines the classical steepest descent method and the two- 
point step size gradient method IZTll . to solve the nonlinear 
hxed-point problem of EM. It only requires the EM updating 
scheme and can be readily implemented as an ’off-the-shelf’ 
accelerator. Since MM is a generalization of EM and the 
update rule is just a hxed-point iteration, SQUAREM can be 
readily applied to MM algorithms. 

Let Fmisl(’) denote the nonlinear hxed-point iteration map 
of the MISL algorithm: 

=Fmisl(x('=)), (52) 

whose form can be dehned by the following equation: 

x(fe+l) ^ g7arg(-A(Diag(p('“))-pW,,I-JV=l)A"x(''))^ ^33^ 

where and pmlx are the same as in Algorithm |2] and 
functions e*-'^ and arg(-) are applied element-wise to the 
vectors. With this, the steps of the accelerated-MISL based 
on SQUAREM are summarized in Algorithm [3 

A problem of the general SQUAREM is that it may violate 
the nonlinear constraints, so in Algorithm |3] the function 
g 7 arg( ) been applied to project wayward points back 
to the feasible region. A second problem of SQUAREM is 
that it can violate the descent property of the original MM 
algorithm. To ensure the descent property, a strategy based on 
backtracking is adopted, which repeatedly halves the distance 
between a and —1 (i.e., a ^ (a — l)/2) until the descent 
property is maintained. To see why this works, we first note 
that X = g 7 arg(x('')- 2 ar+a^v) = X 2 if « = — 1. In addition, 
since ISL(x 2 ) < ISLjx^^^) due to the descent property of 
original MM steps, ISL(x) < ISL(x('')) is guaranteed to hold 
as a —— 1. In practice, usually only a few backtracking steps 
are needed to maintain the monotonicity of the algorithm. 


B. Acceleration via Backtracking 


As mentioned earlier, the slow convergence of MISL could 
be due to the double majorization, which may have resulted 
in a majorization function that is a very loose approximation 
of the original objective function. Thus to accelerate the 
convergence, one possibility is to hnd a better approximation 
of the objective at each iteration. Notice that to preserve 
the descent property (|25ll of the majorization-minimization 
method, it only requires m(x, > /(x) at x = x(^+^\ 
i.e., at the minimizer of rt(x, x^^^), rather than requiring 
rt(x, x^*)) to be a global upper bound of /(x). With this in 
mind, we consider the following approximation function of 


/W =Ep=i at xW: 


Algorithm 3 Accelerated-MISL. 
Require: sequence length N 
1: Set k = 0, initialize x’^°\ 

2: repeat 

3: Xi = FmiSl(x*'^^) 

4: X2=FmiSl(xi) 

5: r = Xi - X^'") 

6 : V = X 2 — Xi — r 

7: Compute the step-length a = 

g. X = giarg(xf*''-2Qr+Q^v) 

9: while ISL(x) > ISL(x('')) do 

10: a {a — 1)/2 

11- X = g7arg(x(^'-2ar+a^v) 

12: end while 

13: x('=+^) = X 

14: fc ^ fc -f 1 

15: untU convergence 


ul(x,x('=)) =4Re (x^A (Diag(p('=)) - Ll) A^x^''^) 

2N 

+ 8N^L - 3^ 

(54) 

where L > 0 needs to be chosen such that rti(x, x^^^) > /(x) 
at the minimizer of Mi(x, x^^^) over the constraint set C = 
{x G C^l \xn\ = l,n = 1,... ,N}. It is easy to see that the 
minimizer of ul (x, x^^^ ) over C is given by 

X* = giarg(A(LI-Diag(p('“)))A"x(''))^ ^33^ 

where functions and arg(-) are applied element-wise to 
the vectors. Thus we need to choose L such that 

f{xl). (56) 

It is worth noting that by taking L = pmlx + N'^, the function 
wz,(x,x^^^) is just the majorization function of /(x) in ( |39] |. 
Hence we know that the condition (l56l l can be guaranteed if 
L is greater than pmlx + N'^- To hnd a better approximation, 
here we adopt a backtracking strategy to choose L. More 
specihcally, at iteration k we choose L to be the smallest 
element in {pmax + (2*'“ — l)A^}ifc=o,i,... such that condition 
(l56l l is satished. Note that by choosing L in this way, the 
resulting function ul (x, x^^^ ) is not guaranteed to be a global 
upper bound anymore, but the monotonic property of the 
algorithm is ensured. We call this algorithm as backtracking- 
MISL and the steps of the algorithm are summarized in 
Algorithm |4] 

V. ISL Minimization with spectral constraints 

In some applications, for example in cognitive radar ll28l . 
apart from designing sequences with good correlation prop¬ 
erties they need to satisfy some spectral constraints like the 
spectral power in some frequency bands should be lower than 
some specihed level. In practice, designing such sequences is 
a challenge, there are few methods available in the literature 


afx('=) 







Algorithm 4 Backtracking-MISL. 

Require: sequence length N 

1 

Set k = 0, initialize x^^k 

2 

repeat 

3 

p('=) = 

4 

Pmlx = maxp{p^^^ : p = 1,..., 2A} 

5 

Set ife = 0 

6 

repeat 

7 

L = pL^lx + (2*'^ - 1)A 

8 

X* = gjarg(A(LI-Diag(p(*^)))A''fx('=)) 

9 

ife *fc + 1 

10 

until ml(xJ^,x('=)) >/(x^) 

11 

^{k+l) _ X* 

12 

k k + 1 

13 

until convergence 


dealing with this problem, see ||20|. In this section, we 
will show how MISL can be adapted to design sequences 
with spectral constraints and we call the resulting algorithm 
spectral-MISL. 

The spectral constraints such as the power in some band, 
denoted hereafter by the set of indices C [1,27V], should 
be lower than a threshold can be expressed as: 

|af xp < e, (57) 

fcen 

where e denotes some pre-specified threshold. Then the prob¬ 
lem of minimizing ISL subject to some spectral constraints 
can be expressed as 

N-l 

minimize 

fe=i (58) 

subject to |a:„| = 1 , n = 1 ,..., 

For a given e > 0, we can always find a A such that the 
problem (l58l l can be transformed into the following equivalent 
problem: 

N-l 

minimize ^ + A^ jaf xp 

fc=i ken ^ ^ 

subject to |a:„| = 1 , n = 1 ,..., 

As shown before in (| 8 ]l, the problem can be expressed as 

2N 

minimize ^ [af xx^a^ - A^] -f A ^ jaf xp 

^ ken ^ ^ 

subject to |a;„| = 1 , n = 1 ,..., A^. 


From here on we can follow the derivation of MISL to derive 
the spectral-MISL algorithm for problem (l60l l. For instance, 
after the first majorization at a given point x*^^\ the majorized 
problem is 


2N 

minimize 2^^|a^x*^^)| x^apa^x 

p=i 

_4Ar2xt^xW(xW)^x 
subject to |x„| = 1, n = 1,..., A^, 


ken 


( 61 ) 


which can be rewritten as 

2N 

minimize 2 ^ x^ ap a^ x — lAf^x^x^^^ (x^^^J^x 

p=i 

subject to |x„| = 1, n = 1,..., A^, 
where 

|2 

4- A /‘.t n 1= W. 

(63) 


(62) 


f,(k) ^ J+ V2, p e 12 

otherwise. 


Pp ' = 


ja^xWp 


By defining = 
we get the problem 


■_(fe) _(fe) 

Pi \ ---,P2N 


1 T 


and A = [ai,... ,a 2 Ar], 


minimize x^ ^ADiag(p(^^)A^ — 2A^x(^^ 
subject to |a;„| = 1, n = 1,..., A, 



which has the same form as the problem (l3^ . So we can 
follow the same steps as before and derive the spectral-MISL 
algorithm, the main steps of which are listed in Algorithm 
|5] The acceleration schemes described in section |IV] can be 
readily applied to spectral-MISL as well, however we will not 
elaborate on that here. 


Algorithm 5 Spectral-MISL. 

Require: sequence length N, index set 12 and A. 

1 : Set k = 0, initialize x^^^. 

2: repeat 

3 . pW ^ + V 2 , pef 2 

1 , otherwise 

4: Pmlx = maxplp^^'^ : p = 1,..., 2A} 

5: y = - A (Diag(p('=)) - pLIxI - A^^l) A^x^'^) 

6: n = 1,...,A 

7: k k + 1 

8 : until convergence 


VI. Numerical Experiments 

To compare the performance of the proposed MISL algo¬ 
rithms with existing ones and to show the potential of MISL 
to design sequences for various scenarios, we present some 
experimental results in this section. All experiments were 
performed on a PC with a 3.20GHz i5-3470 CPU and 8 GB 
RAM. 

A. ISL Minimization 

In this subsection, we present numerical results on apply¬ 
ing the proposed algorithms to design unimodular sequences 
with low ISL and compare the performance with the CAN 
algorithm El, which is known to be computationally very 
efficient. The Matlab code of CAN was downloaded from the 
websit43 of the book 141 . 

'http://www.sal.ufl.edu/book/ 
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Figure 1. Merit factor (MF) versus sequence length. Each curve is averaged 
over 100 random trials. 



Figure 2. Average running time versus sequence length. Each curve is 
averaged over 100 random trials. 


We first compare the quality, measured by the merit factor 
(MF) debited in (|2|i (recall that the larger the better), of the out¬ 
put sequences generated by different algorithms. In this exper¬ 
iment, for all algorithms, the stopping criterion was set to be 
|lSL(x('=+i))-ISL(x('=))|/max(l,ISL(xW)) < and 
the initial sequence was chosen to be 

where {On]n=i ^re independent random variables uniformly 
distributed in [0,1]. Each algorithm was repeated 100 times for 
each of the following lengths: N = 2®, 2®,..., 2^®. The aver¬ 
age merit factor of the output sequences and the average run¬ 
ning time are shown in Figs. [Hand |2] for different algorithms. 
From Fig. [T] we can see that the proposed backtracking-MISL 
and accelerated-MISL algorithms can generate sequences with 
consistently larger MF than the CAN algorithm for all lengths; 
and at the same time they are faster than the CAN algorithm 
as can be seen from Fig. [2] especially the accelerated-MISL 
algorithm. Taking both the merit factor and the computational 
complexity into consideration, the accelerated-MISL seems to 
be the best among the three algorithms. The correlation levels 
of two example sequences of length V = 512 and N = 4096 
generated by the accelerated-MISL algorithm are shown in 
Fig. a where the correlation level is debned as 


correlation level = 20 logj^g 


T’O 


, fc = 1 - iV,..., TV - 1. 


Next, we perform an experiment to show the different 
convergence properties of the proposed accelerated-MISL al¬ 
gorithm and the CAN algorithm. For both algorithms, we 
brst initialize them with the same randomly generated se¬ 
quence. Then we initialize CAN with the output sequence 
of accelerated-MISL and initialize accelerated-MISL with the 
output sequence of CAN, and run the two algorithms again. 
The evolution curves of the ISL in two different random trails 
with N = 32 are shown in Fig.|4|a) andHJb). From the bgures, 
we can see that when initialized with the output sequence of 
CAN, the accelerated-MISL can further decrease the ISL a 
little bit, which means the point CAN converges to is probably 
not a local minimum. While in contrast, when initialized with 


N = 512 



k 


Figure 3. Correlation level of two sequences of length N = 512 and 
N = 4096, generated by the accelerated-MISL algorithm. 


the output sequence of accelerated-MISL, CAN increases the 
ISL, which means the point accelerated-MISL converges to 
is probably a local optimal and also tells the fact that CAN 
does not share the monotonicity of MISL. In addition, we can 
also see that when initialized with different sequences, CAN 
and accelerated-MISL may converge to different points. Thus, 
initializing the algorithms with good sequences, e.g. Frank 
sequence na, Golomb sequence US, can probably improve 
the performance of the algorithms. 

Finally, we present an example of applying the accelerated- 
MISL algorithm to design unimodular sequences with 
impulse-like periodic autocorrelations. Recall that in the pe¬ 
riodic case, unimodular sequences with exact impulsive au¬ 
tocorrelation (i.e., zero ISL) exist for any length N. The 
periodic correlation levels of two generated sequences of 
length N = 256 and N = 1024 are shown in Fig.jH where the 
periodic correlation level is debned as in the aperiodic case. 
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Figure 4. Evolution of the Integrated sidelobe level (ISL) in two different 
random trails with TV = 32. 



k 

Figure 5. Periodic correlation level of two sequences of length TV = 256 
and TV = 1024. 


From Fig. 0 we can see that the periodic correlation levels 
of the two sequences are indeed very close to zero (lower than 
-200dB) at lags fc 7 ^ 0. 

B. Spectral-MISL 

In this subsection, we consider the design of a sequence 
of length N — 1000 with low autocorrelation sidelobes and 
at the same time with low spectral power in the frequency 
bands [f, f) U [^,7r) U To achieve this, we set 

the index set in spectral-MISL as 17 = {250,..., 499} U 
(750,..., 999} U (1500,..., 1749} accordingly and tune the 
parameter A to adjust the tradeoff between minimizing the 
correlation level and restricting the power in the pre-specified 
frequency bands. Fig. |6] shows the spectral power of the 
output sequence generated by spectral-MISL when initialized 
with a random sequence and with parameter A = 10^. The 
correlation level of the same output sequence is shown in 
Fig. I 2 ] We can see from the figures that the power in the 
pre-specified frequency bands has been suppressed and the 
correlation sidelobes are still relatively low. 



_4Q I-U- MM-U- U-LI 1-1—I- 

0 71/2 n 3 k/1 271 

angular frequency (rad/s) 


Figure 6. Spectral power of the designed sequence of length TV = 1000, 
generated by spectral-MISL with A = lO'^. 

VII. Conclusion 

We have presented an efficient algorithm named MISL 
for the minimization of the ISL metric of unimodular se¬ 
quences. The MISL algorithm is derived based on the general 
majorization-minimization scheme and the convergence to a 
stationary point is guaranteed. Acceleration schemes that can 
be used to speed up MISL have also been considered. Numer¬ 
ical results show that in the case of aperiodic autocorrelations 
the proposed MISL algorithm can generate sequences with 
larger merit factor than the state-of-the-art method and in the 
case of periodic autocorrelations, MISL can produce sequences 
with virtually zero autocorrelation sidelobes. A MISL-like 
algorithm, i.e., spectral-MISL, has also been developed to 
minimize the ISL metric when there are additional spectral 
constraints. Numerical experiments show that the spectral- 
MISL algorithm can be used to design sequences with spectral 
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Figure 7. Correlation level of the designed sequence of length N = 1000, 
generated by spectral-MISL with A = lO"^. 

power suppressed in arbitrary frequency bands and at the 
same time with low autocorrelation sidelobes. All the proposed 
algorithms can be implemented by means of the FFT and thus 
are computationally very efficient in practice. 

Appendix 

A. Proof of X^U^)=2N^ 

Proof: We first recall that 

2N 

= y]vec(Ap)vec(Ap)'^, (65) 

p=i 

where Ap = apa^, ap = [1, ^ 

and Wp = - 1), p = The 

(m, n)-th element of matrix Ap is given by 

which is also the (m + (n — 1) A)-th element of 
vec(Ap). Then the (mi + (ni — 1)A, m 2 + (n 2 — 1) A)- 
th element of the matrix vec(Ap)vec(Ap)^ is given 
by [vec(Ap)vec(Ap)^]^^^(^^_^^^_^^^(^^_^^^ 
gja;p(mi-ni-m2+n2)^ = 1,..., A, Tli = 1,...,A, m 2 = 

1.. .., A, 712 = 1,..., A. Then it is easy to see that 

2N 

rasl _ „j<^p(mi-rai-m2+"2) 

l^Jmi+(ni-l)Af,m2 + (n2-l)W / y ^ 

p=l 

2A, mi — 1712 =711—772 
0, otherwise. 

(66) 

We can see that $ is a real matrix. 

To prove Aniax(^) = 2A^, we will show that 

(2A^I — $) X > 0, Vx G and the equality is 

achievable for some x ^ 0. Let us define Iq = {1,A + 

2. . .., 7v2 - A - 1, A^}, I_i = {2, A + 3,. .., (A - 1)A}, 
Ii = {A + 1,2A + 2,...,A2-1 },...,Ii_iv = {A} and 


^N-l = {A^ — A + 1}. It can be shown that 

x'T (2A2I - $) X 

( N-l N-1 \ 

H H (xi-xjf+ ■ 

fc—1—N k—l — NiGXk J 

(67) 

From (l67T i. it is easy to see that x^ (2A^I — $) x > 0, Vx G 

1\T^ 7V7^ 

R and the equality holds for all x G R with the 
following structure; 

Xi Xj , J ^ Tq 

Xi = 0, otherwise. 

The proof is complete. ■ 
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